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Abstract 

When analyzed in terms of the Symanzik expansion, lattice corre¬ 
lators of multi-local (gauge-invariant) operators with non-trivial con¬ 
tinuum limit exhibit in maximally twisted lattice QCD “infrared di¬ 
vergent” cutoff effects of the type a?^/{rn^)^, 2k > h > 1 {k,h inte¬ 
gers), which tend to become numerically large when the pion mass gets 
small. We prove that, if the action is 0(a) improved d la Symanzik 
or, alternatively, the critical mass counter-term is chosen in some “op¬ 
timal” way, these lattice artifacts are reduced to terms that are at 
worst of the order {a?/, k>l. This implies that the contin¬ 
uum extrapolation of lattice results is smooth at least down to values 
of the quark mass, mq, satisfying the order of magnitude inequality 
mq > a^A^cD- 


1 


1 Introduction and main results 


Although in lattice QCD at maximal twist (Mtm-LQCD) 0(a) discretization 
effects (actually all 0(a^*'’''^), fc > 0 effects) are absent or easily eliminated [H 
laEi, it turns out that, because of the explicit breaking of parity and iso¬ 
spin symmetry present in the action, correlators are affected by dangerous 
artifacts of relative order a'^^, k > 1, which are enhanced by inverse powers 
of the (squared) pion mass, as the latter becomes small. Terms of this kind 
may have been seen in the simulations reported in ref. 

When analyzed in terms of the Symanzik expansion, the expectation val¬ 
ues of multi-local (gauge-invariant) operators, exhibit, as 0, 

what we will call “infrared (IR) divergent” cutoff effects with a behaviour of 
the form ^ 


{O) 


L 

mq 


{O) 


. ^2k 

r., ^ ® 

mq 1 ^(m^)^ 



2k > h > 1 {k,h integers), 


( 1 . 1 ) 


where we have assumed that the lattice correlator has a non-trivial continuum 
limit. 

It is important to remark that the existence of pion poles in the Symanzik 
expansion of a lattice correlator in no way means that the latter diverges as 
rriq —>■ 0, but only that the Symanzik expansion may break down if the order 
of magnitude inequality > a is violated |5]. 

We will prove that, if the action is 0(a) improved d la Symanzik by in¬ 
troducing the lattice clover term ^ (with its non-perturbatively determined 
csw coefficient [Zj), or the critical mass counter-term is chosen in some “op¬ 
timal” way, these artifacts will be reduced to terms that are at worst of order 
a^(a^/m^)^“^, k >1. 

The idea that a suitable definition of critical mass exists which can lead 
to a smoothing out of chirally enhanced lattice artifacts or perhaps be of 
help in getting improvement was already put forward in the context of chiral 
perturbation theory (xPT) in refs. jH] and [Hj, respectively. 

A crucial consequence of the analysis we present in this paper is that the 
strong (order of magnitude) inequality 


iriq > uAqqj^ , (1.2) 

invoked in ref. j2] in order to have the phase of the chiral vacuum driven by 
the quark mass term and not by the (twisted) Wilson term, can be relaxed 
to the weaker relation 

rriq > a Aq^^ , (1-3) 

^ Often for short we do not include powers of Aqcd required to match physical dimen¬ 
sions. 
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before large cutoff effects are possibly met when the quark mass is lowered 
at hxed a. The bound (Q is fairly weak as it allows simulations in a region 
of quark masses that correspond to rather light pions (with masses close to 
around 200 MeV for typical present-day lattice spacings). 

In order to make this observation useful in practice, it is however neces¬ 
sary to have a handle on the numerical coefficients entering the bound m- 
We will argue that one can monitor the “safe” region of quark masses by 
measuring the mass of the charged pion as a function of rUg. It can be 
shown, in fact, that, once the leading “IR divergent” cutoff effects have been 
eliminated (and under the assumption that chiral symmetry is spontaneously 
broken in the continuum), the lattice squared mass of the charged pion, for 
sufficiently small rriq, is a linear function of (with only small O(a^) and 
0(mq) distortions) at least until the pion mass remains within the bound 
hinted at by eq. (USD- 

The situation one finds in this kinematic regime is thus reminiscent of 
the ideal continuum case, where the Gell-Mann-Oakes-Renner (GMOR) re¬ 
lation HDl ensures the proportionality of to the mass of the quark. Armed 
with the lattice analog of the GMOR relation one can exploit the previous 
considerations to argue that, vice-versa, at fixed a, rUg should not be de¬ 
creased beyond the point where non-linearities show up. 

The effectiveness of Mtm-LQGD in removing 0(a) discretization errors 
was successfully tested in quenched simulations in refs. HH. The ability 
of the optimal choice of the critical mass in diminishing the magnitude of 
lattice artifacts at small quark mass has been beautifully demonstrated in 
the recent works of refs. [nmnnn- Results of unquenched simulations na, 
though encouraging, are still rather preliminary due to the complicated phase 
structure cni of Wilson fermions when mq/AqcD is numerically comparable 
to (or smaller than) a^Aqqj^. 

The plan of the paper is as follows. In Section [2l we analyze the form of 
the Symanzik expansion of lattice correlators beyond 0(a) and explain why 
and how “IR divergent” cutoff effects arise in this context. In Section IHl we 
discuss two ways of removing all the leading “IR divergent” cutoff effects and 
we describe the structure of the left-over “IR divergent” terms. In Section |3] 
we derive the lattice GMOR relation and we discuss how it can be used to 
monitor the lowest value of the quark mass that can be safely employed if the 
lattice spacing is held fixed. Finally in Section El we collect some remarks on 
the peculiar structure of the lattice artifacts affecting lattice hadron energies 
and the pion decay constant. Gonclusions can be found in Section El In 
an Appendix we prove that automatic 0(a^^’''^) improvement of the lattice 
expectation values of parity-even operators holds in Mtm-LQGD, by using 
an argument which does not rely on the spurionic transformation r —^ —r. 
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hence on the r-parity properties of the critical mass. 


2 Symanzik analysis of “IR divergent” cutoff artifacts 


The expression of fermionic action of Mtm-LQCD is given in the physical 
qnark basis 121 ini by the formnla 


oL 

^tmF 


= a 




7-V + mg - ^75X3 


( 


r 

a-V* 


■V + M,^,(r 


r 


(X 


( 2 . 1 ) 


where is a lattice fermion donblet, r is the Wilson parameter and 
is the bare qnark mass, i.e. the real parameter which in the continnnm will 
provide a non-vanishing mass to the pion. is any sensible (from the point 
of view of renormalization theory) “estimate” of the critical mass which, in 
order to match the r-parity properties of the Wilson term, should be taken 
as an odd function of r j2] . 

The study of cutoff artifacts affecting lattice correlators in Mtm-LQCD 
can be elegantly made in the language of the Symanzik expansion. A full 
analysis of cutoff effects beyond 0(a) is of course extremely complicated. 
Fortunately it is not necessary, if we limit the discussion to the terms that 
are most strongly enhanced as the quark mass is decreased. This analysis 
will be carried out in the next section. 

For completeness we show in an Appendix that automatic 0(a) (actually 
0 (a 2 fc+i), > 0) improvement directly follows from the symmetry of the 

lattice theory under the transformation P x x (m^ ^ —^q), where P is 
parity and the transformation x [rriq —>■ —rriq) counts the parity of the 
overall dimension of an operator. The explicit dehnition of the parity opera¬ 
tion P and the transformation Pd can be found in eqs. (EH) and (Toil . The 
proof we give does not involve any change of sign of the Wilson parameter, 
thus it is independent of the r-parity properties of the critical mass unlike 
the argument developed in ref. j2]. 


2.1 The Symanzik LEEA of Mtm-LQCD 

With reference to the fermionic lattice action m, the low energy effective 
action (LEEA) of Mtm-LQCD, S'sym, can be conveniently written in the form 


‘S'svm 


^We wish to thank M. Liischer for a stimulating discussion on this issue during the 
“Twisted Mass Lattice Fermions” workshop held in Villa Mondragone (Frascati - Italy) 
on March 14-15, 2005. 


3 









where £4 = ^^triF-F) + 'ipH-D + is the target continuum QCD 

^9q ^ 

Lagrangian density. A number of interesting properties enjoyed by the above 
LEEA can be proved which we summarize below. 

• Lagrangian density terms of even dimension, i 2 k, in eq. (I2.2|l are parity- 
even, while terms of odd dimension, i 2 k+i, are parity-odd and twisted in 
iso-spin space. Thus the latter have the quantum numbers of the neutral 
pion. These parity properties follow from dimensional arguments and the 
invariance of the LEEA of Mtm-LQCD under the transformation P x x 
{niq —>■ —mq), inherited from the correspondent invariance of the lattice 
theory (see Appendix). 

• The term of order a in eq. (El, 4, is given (after use of the equations 
of motion of continuum QCD) by the linear combination 

4 = S5,SW ^5,SW + <^5,^2 + 55,e 4,e , (2.3) 

X — — — 

■ ^]*75'r3^ , 4,m2 = = A^cd^^sDV' , (2.4) 

where the coefficients S^^sw, 4 ,m 2 and 4 ,e are 0 ( 1 ) dimensionless quantities, 
odd in r. The operator £5 5 arises from the need to describe order a uncer¬ 
tainties entering any non-perturbative determination of the critical mass and 
goes together with £ 5 ^sw- Both i^^sw and 4,e could be made to disappear 
from eq. (El by introducing in the Mtm-LQCD action the SW (clover)-term 
with the appropriate non-perturbatively determined csw coefficient and at 
the same time setting the critical mass to its correspondingly 0 (a) improved 
value. 

• Higher order ambiguities {k > 1) in the critical mass will be described 
by terms proportional to odd powers of a, more precisely of the kind 

4+2fc-|-l,e 4+2A:-|-l,e = 4+2A:+l,e (Aqcd)^^"''^ • (2.5) 

The structure of these terms, which will all contribute to £odd (see eq. (El 
below), follows again from obvious dimensional arguments plus the invariance 
of the lattice action under the transformation P x x {'mq —'mq). 

2.2 Describing Mtm-LQCD correlators beyond 0(a) 

We are interested in the Symanzik description of the lattice artifacts af¬ 
fecting connected expectation values of n-point, multi-local, multiplicative 
renormalizable (m.r.) and gauge-invariant operators 

n 

0{Xi,X2,...,Xn) = Y[0j{Xj)=0{x), Xi ^ X2 ^ ^ Xn, ( 2 . 6 ) 

J=1 
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which we take to have continuum vacuum quantum numbers, so as to yield 
a result that does not trivially vanish as a —>■ 0. In particular, in order 
to ensure automatic 0(a) improvement j2j we will assume that O is parity 
invariant in which case its Symanzik expansion will contain only even powers 
of a. Schematically we write 


(0(a;)) = ([0(a;) + AoddO(x) + , (2.7) 

ruq rriq 

where for short we have introduced the compact notations 

( 2 . 8 ) 


-^odd — X! , 

k=0 


-^even E 

k=l 


The operators AoddO (resp. AevenO) have their origin in the need of regu¬ 
larizing terms where a parity-odd (resp. parity-even) product of £odd and/or 
£even insertions comes in contact with some of the points where the local 
operator factors appearing in O are concentrated. AoddO (resp. AgvenO) 
counter-terms have an expansion in odd (resp. even) powers of a. We recall 
that they can be viewed as the n-point operators necessary for the on-shell 
improvement of O nam 

For the purpose of this discussion we imagine expanding the r.h.s. of 
eq. (EH) in powers of Jd^yCoddiv) and/or £even(2/)- Terms with j and/or 
j' insertions of the hrst and/or the second of these factors will generate in 
the Symanzik expansion h-fold 1/m^ pion poles with 1 < h < j + j'. 


2.3 Pion poles and “IR divergent” cntoff effects 

Although a complete analysis of all the “IR divergent” cutoff effects is very 
complicated, the structure of the leading ones {h = 2k in eq. (jl.ljl l is 
rather simple, as they only come from continuum correlators where 2k factors 
/ d'^y Codd{y) are inserted. More precisely the leading “IR divergent” cutoff 
effects are identihed on the basis of the following 

Result: In the Symanzik expansion of {0{x))\^^ at order {k > 1) 
there appear terms with a 2/c-fold pion pole and residues proportional to 
|(r2|£odd|7r°(0))p*', where (r2| and |7r°(0)) denote the vacuum and the one- 

state at zero three-momentum, respectively. Putting different factors 
together, each one of these terms can be seen to be schematically of the form 
(recall £odd = 0(a)) 


/I \ 2k 

(/ft) (^-("^-?))^^-^[O;{vr°(0)}2fc] 


cont 

rriq 


(2.9) 
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where we have generically denoted by A^[ 0 ; { 7 r‘^( 0 )} 2 fc] the 2 /c-particle matrix 
elements of the operator O, with each of the 2k particles being a nentral pion 
at zero three-momentum and we have introduced the short-hand notation 




(f2|£odd|vr°(0)) 


cont 

rriq 


( 2 . 10 ) 


The proof of the above Result can be obtained on the basis of 1) the 
general theorems of quantum held theory governing the appearance of poles 
in correlators ^HI, 2 ) the notion of “partially disconnected” diagrams which 
come about when reducing pairs of pions from in and out states of a multi- 
pion matrix element, 3) the observation that the (2A;-|-l)-pion matrix elements 
of £odd can all be shown to be proportional to in the chiral limit, through 
a repeated use of soft pion theorems HHI, in which pions are successively 
reduced out. 

Less “IR divergent” cutoff effects (those with h strictly smaller than 2k 
in eq. fjl.lji l come either from terms with insertions of / (RyCeveniv) and/or 
pairs of / (RyCo^^iy) factors, or from contributions of more complicated in¬ 
termediate states other than straight zero three-momentum pions or from 
both. In all cases one gets terms with extra a? powers not “accompanied” 
by an equal number of factors. 


3 Reducing “IR divergent” cutoff artifacts 

We have shown in the previous sections that close to the chiral limit the 
most “IR divergent” discretization artifacts affecting Mtm-LQCD correla¬ 
tors at order are proportional to 2k powers of (see eq. (I23D). Since 
>^^odd = 0 ^ 5 - 1 - 0{R), this also means that at leading order in a each multiple 
pion pole residue is proportional to |(r 2 |£ 5 | 7 r°( 0 ))p^. It is an immediate con¬ 
clusion of this analysis that all these dangerous cutoff effects can be removed 
from lattice data if we can either eliminate £5 from the Symanzik LEEA of 
Mtm-LQCD or set to zero. Actually, in the last case a somewhat weaker 
condition is sufficient. Indeed, we will see that it is enough to reduce to a 
quantity of order am^. 


3.1 Improving the Mtm-LQCD action by the SW-term 

The obvious, field-theoretical way to eliminate £5 from the Symanzik LEEA 
of Mtm-LQCD consists in making use of the 0(a) improved action 


slmF = + ^<1 + 
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+ Mi;(r) + jCsw{r)[a-F]^)]^\x ), (3.1) 


i 

4^ 


where csw is fixed in the appropriate non-perturbative way [7j and is 
an improved estimate of the critical mass. 

In this situation the lattice correlation functions of the theory will admit a 
Symanzik description in terms of a LEEA where the operators i^^sw and 
are absent, and £5 will be simply given by £5,^2 (see eq. Q). The left-over 
contributions arising from the insertions of in will yield terms 

that are at most of order [arrig/rn^y^ ~ (am^)^^, hence negligible in the 
chiral limit. It is instead the next odd operator in the Symanzik expansion, 
which comes into play. 

A detailed combinatoric analysis based on the structure of the non-leading 
“IR divergent” cutoff effects reveals that the worst lattice artifacts left behind 
in correlators after the clover cure are of the kind a^(a^/m^)^“^, k > 1. 


3.2 Optimal choice of the critical mass 


The alternative strategy to eliminate all the leading “IR divergent” cutoff 
effects consists in leaving the Mtm-LQCD action of the form but fixing 
the critical mass through the condition 


lim 

mg —*■0 + 


lim 

mq—>0+ 


(f]|£oddk°(0)) 


cont 

niq 


0. 


(3.2) 


The meaning of this condition is rather simple. It amounts to fix, for 
each k > 0, the order contribution in the critical mass counter-term, 

—M°[^ ^ in such a way that its vacuum to one-pion state matrix el¬ 

ement compensates, in the limit of vanishing quark mass, the similar matrix 
element of the sum of all the other operators making up £4+2fc+i- 

In the next section we present concrete procedures designed to implement 
the condition (IS2D in actual simulations. To avoid confusion with the values 
the lattice action parameters will take in the successive steps of a simulation, 
we will provisionally put a bar over the symbols representing the values of the 
quark and the corresponding pion state mass while we develop the argument 
for the “optimal” determination of the critical mass. 


3.2.1 Lattice estimate 

We want to show how eq. dH, which has to do with a matrix element 
defined in the continuum theory, can be translated into a lattice condition. 
To this end let us consider the lattice correlator 




Xci 7^ 0, 


(3.3) 
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where Vq = is the vector current with iso-spin index 2 and = 

'075^'*/' is the pseudo-scalar quark density with iso-spin index 1. In the 
continuum the correlator (HI owing to parity vanishes for any value of rUg, 
and we have 


{QvP\0)) 


cont 

friq 


0 , 


(3.4) 


where Qy = fdxVQ(x,t) is the iso-spin 2 vector charge. On the lattice the 
breaking of parity (and iso-spin) due to the presence of the twisted Wil¬ 
son term makes the correlator (Q non-vanishing by pure discretization 
effects. Extending to parity violating correlators the arguments developed in 
Sect. 12.11 one gets for its Symanzik expansion 


(Q(a:,0)) _ = 

^' ran 

X ^ 

= J dx{(AoddQ(2:, 0)) 


(Q(a:,0) Jd'^y Coddiy)) 


cont 

fhq 


+ ..., (3.5) 


where dots represent terms with higher order insertions of / d'^y Codd{y) and/or 
Jd'^y Ceveniy) and AoddQ has the expression 


AoddQ(a:,0) = 


(3.6) 


= a 


Tl 


h^o^/’75y^(a^)^75^V’(0) -Fr7m,V’7o75yV'(3^) V’75yV’(0)J + 0{a^'' 


In eq. dnn) rj and r] are appropriate dimensionless coefficients, odd in r. We 
recall that in the language of the Symanzik improvement program the term 
we have explicitly written down in eq. dSH) is the standard operator necessary 
for the on-shell 0(a) improvement of Q(x, 0). In writing the expansion ()3.5|1 
use has been made of the continuum relation dH- 

For the rest of the argument it is important to remark that from higher or¬ 
der insertions of ^d'^y Coddiy) arbitrarily high powers of the ratio ^T,ifhq)/ffi^ 
will be generated in the r.h.s. of eq. (13.5|) 

At large times [t S> 1/Am, where Am is the difference between the mass 
of the hrst excited state with the pion quantum numbers and that of the 
pion), we may write more explicitly the r.h.s. of eq. (j3.5jl in the form 


re 


lim 

t^+oo L 2m^ 


-1 




amqT] (0|Ao|7r (0 ))C'ssl + 


-am^r]{VL\P^\T^\Q))CsL - ^m^CL ] (7r^(0)|P^|O) 

2m£ J ^<1 


(3.7) 


with all corrections of order or higher encoded in the coefficient factors 


Cx — 1 + 


A) 


OO 


/ ^7r('^q) V / ^ ^ 
V ^2 J \ffi2j 


X = L,SL,SSL. (3.8) 


j'+.^=even ’ 


















Several observations about the r.h.s. of eq. (EH) are in order here. 1) Since 
mq(r2|y4Q|7r^(0))|5^“* is a quantity of order the hrst term is completely 
immaterial to the present analysis. 2) The factor rh-„^ in front of the sec¬ 
ond term comes from the time derivative in eq. (jd.tijl . 3) The factor 2rn^ 
in the third term arises from the chain of relations (7r°(0)|(5y|7r^(q))|^’^* = 
(7r°(0)|7r°(q))|5^“* = (27r)^2m^|5^“M^(q). 4 ) Concerning the factors Cl, Csl 
and CssL, they differ from unit by three different types of lattice artifacts, 
i) The terms with /c = £ = 0 in Cl represent the the leading “IR divergent” 
corrections contributing to the correlator Sx(2(^) 0))lm • ii) Sub-leading 
“IR divergent” corrections are generated by insertions of / d^yCe^eniy) and/or 
pairs of / d'^yCoAdiy) factors. They contribute extra /m^ powers, iii) Fi¬ 
nally there are “IR hnite” corrections stemming, among others, from contri¬ 
butions of intermediate states other than zero three-momentum pions. 

With all these premises, in order to fix the critical mass so as to have in 
the continuum (see eq. dSH and the dehnition (ITTUld 


lim ^^(mg)=0, 

mg —*■0 + 


(3.9) 


one may think of proceeding on the lattice in the following way. 

Given a first estimate of the critical mass, say M/,., consider the lattice 
action dj where the critical mass has been momentarily replaced by the 
expression M/,. -t- 5fh. In order to implement eq. (uni) we must compute the 
lattice quantity in the l.h.s. of eq. dSH and identify the optimal value of the 
critical mass, M/f*, as the limiting value of -1- 5m at which the 0(a°) 
quantity (see eqs. dSH and dSH) 


Airriq, M/j, + 5m] t) 


[m 


(Q(x,0)) 

a ^ 




t > 


1 

Am 


(3.10) 


vanishes as ffiq is extrapolated to smaller and smaller values from the region 
where ffiq > a. 

Numerically there can be various ways to do this. One possible strategy 
is to start from a value, frSq \ of the quark mass such that the order of 
magnitude inequality > a holds. A hrst determination, 5m^^\ of 5m 

can be obtained by enforcing at large t the condition 

+ 5m^^^]t) = 0 ^ , (3-11) 

^ Am 

which in turn yields a hrst estimate, of the matrix element in ()3.9|) . Solv¬ 
ing iteratively the non-linear equation dHH), one can write from eqs. dSH 
and ()3.8j) 

= -a(rn«)"r,(ll|/’V'(0)) + , (3.12) 
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where the hrst term is the solution of the linearized form of (ITTTD obtained 
by setting to unit the factors Cl, Csl and Cssl in eq. (EUD- The second 
term, is the correction due to all the other terms and, in particular, to 

the higher powers of To estimate the magnitude of we plug back the 
ansatz (imi in (jSIH). The structure of the result is 


(^TTT-tt ) 


k=0 


(3.13) 


where the leading O(a^) corrections can be seen to arise from the terms with 
k = 1, j = i = 0 in eq. (EHl) ^ The important point about eq. is that 

can be considered small under a condition, > a^, which is weaker 
than the one > a) we have been using in establishing this result. 

At this point one continues by lowering the quark mass to < fh^g\ 
seeking for the new value, 6rh^^\ of the mass shift which makes vanishing. 
The search must be performed in the neighborhoods of to be sure 

that one remains in a region of the + hm)-plane, where is 

small, so that higher powers of the ratio are even smaller. Proceeding 

in this way the convergence of the Symanzik expansion (HID is not put in 
danger as rfig is decreased. Rather, the initial convergence bound frig > a is 
progressively weakened towards rfig > a^, as signaled by eq. dsm. 

A sequence, of mass shifts is thus determined. If desired, these 

values can be numerically extrapolated to rhg —> 0. The limiting value, 
obtained in this way will allow to identify the “optimal” critical mass 
as the quantity 

. (3.14) 

Notice that in all this procedure no lattice data points are employed where 
the bound fhg > is violated. This caution is necessary because for rhg < 
0 ? large cutoff effects, which are hinted at by uncanceled non-leading “IR 
divergent” terms in the Symanzik expansion, cannot in general be excluded 
and reliable simulations may even be impossible because of metastabilities, 
if e.g. the peculiar lattice phase structure known as the “Sharpe-Singleton 
scenario” HE] is realized. 

The method discussed above may seem unpractical in view of the fact 
that, especially for unquenched simulations, producing data at several values 
of the bare quark mass, as it is necessary to do if one wants to extrapolate 
to fhg —> 0, is computationally rather demanding. We immediately notice, 

^Actually only terms from Cl and Csl are important to this order. Terms coming from 
Cssl are, in fact, negligible because of the extra fhg factor in front of the corresponding 
matrix element. 
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however, that to all practical purposes (see the detailed argument given in 
Sect. Id .2. 21 below 1 one can avoid such an extrapolation and just work with 
the “quasi-optimal” critical mass 

, (3.15) 

where is the smallest of the bare quark masses of interest at the 

lattice spacing one is working and 6m\fnniin is the solution of the equation 
A(m““, t) = 0. 

The idea of estimating the critical mass at the smallest available value of 
rUg was already put forward in the analysis of tm-LQCD performed in yPT 
in refs. iHiini and directly used in the unquenched simulations of ref. na. 
In quenched simulations slightly different versions of the strategy described 
above have been implemented in order to evaluate the (quasi-)optimal value 
of the critical mass. They all turned out to be quite effective in reducing 
cutoff artifacts nacsiiiii- 


3.2.2 Left-over “IR divergent” cutoff effects 

To complete our analysis we have to determine what is the order of magnitude 
of the left-over discretization errors that will affect simulations carried out 
at non-vanishing rOg, when either i) eq. (imi) or ii) eq. (ina for the critical 
mass is inserted in the Mtm-LQCD action m- We will see that in both 
cases the situation will be very much like the one we encountered in Sect. 13.1 1 
where we discussed the case in which the clover term was added to the Mtm- 
LQCD action. 

i) Let us denote by the superscript “ ” lattice quantities computed 

using in the fermionic action the optimal critical mass, (eq. (j3.14jl L The 
formerly leading 2fc-fold pion pole contribution (IZSD now have an expression 
where is systematically replaced by fff^{mg), which represents the 

value of the continuum matrix element |(D|£°j(j|7r°(0))| at the current value, 
nig, of the quark mass, with the parity odd part of the Symanzik LEEA 
when the optimal critical mass is employed. 

Although non-zero at non-vanishing quark mass, we want to show that 
^°^*(mq)/m^ = 0(a). To this end we notice that one can write 


= C'to) + 


dnig 



mg + ... , 


(3.16) 


where from eqs. ()3.12j1 and ()3.13j) one has ^°^*(0) = O(a^). This result 
comes from having extrapolated to frig = 0 from the region rfig > of, 
where = O(a^). Since df^jf^iniq)/dmg\mq=D ~ 0{a), we conclude that 
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is reduced to a mere order am^ quantity, implying, as announced, 

^°P\mg)/ml = 0(a). _ 

ii) If we set instead (eq. ()3.15|l i in the action ()2.1|1 . 

we get (notation should be self-explanatory) 


e 


qopt 


>,)=+ 




drur, 


[ma 




) + ..., (3.17) 


where from eqs. (j3.12jl and ()3.13|) one hnds ^°P*(m™™) = 0(am™“), provided 
> a?. From {nig)/drUglffi^in ~ 0(a) one obtains for any > rh™™ 

the estimate ^l°P^{mq)/m\ = 0(a, am™“/mg) = 0(a). 

The conclusion of this analysis is that both with the optimal as well as 
the quasi-optimal critical mass, the formerly leading “IR divergent” cutoff 
discretization errors are reduced to hnite 0(a^^) contributions. 

This does not mean that all the non-leading “IR divergent” cutoff effects, 
which are of order a^^/(m^)^, 2k > h > 1 {k, h integers), have disappeared 
from correlators. Actually by a non-trivial diagrammatic analysis, based 
on the structure of the non-leading “IR divergent” cutoff contributions, one 
can prove that the most “IR divergent” lattice artifacts left behind after 
using either the optimal or the quasi-optimal dehnition of the critical mass 
(eq. ()3.14j) or eq. (HinD) are reduced down to order a?{a?/m^)^ ^ (fc > 1) 
effects at worse, just like in the case the clover term is employed. Notice, 
however, that in the case where the quasi-optimal critical mass is adopted 
this result holds only for rUg > m™™. 


4 The lattice GMOR relation 

In the previous sections we have shown how the leading “IR divergent” dis¬ 
cretization effects can all be eliminated from the Symanzik expansion of 
lattice correlators. The left-over “IR divergent” terms of the expansion turn 
out to have the structure a? Y.i>i ■ The convergence of this series 

sets the order of magnitude inequality a^ < m\ ~ nig from which one should 
determine the minimal value of the quark mass that (at hxed value of a) can 
be safely simulated before possibly encountering large discretization errors. 

A workable way to numerically estimate this minimal value can be ob¬ 
tained by considering the behaviour of the charged pion mass as a function of 
nig. It turns out, in fact, that in Mtm-LQCD there are Ward-Takahashi iden¬ 
tities (WTI’s) which take exactly the form they have in the formal continuum 
theory. From them a lattice GMOR relation can be derived. 

To see how this works we recall that in Mtm-LQCD the 1-point split axial 
currents, with iso-spin index 6 = 1,2 are exactly conserved in the chiral 
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limit rUg = 0 Hi. This implies the validity of the lattice WTI’s 


( d;A^{x)-2mgPHx) P^m ^ = {SAO)) 


t± _ ;il a2 d± _ J.L^. p =*= ^^2 ,,L 


A^ = Al± lAi , = ^^75 


^x,0 ) 


L^l.L 




II fl " II ■’ - r lu 2 

After integration over space-time, one gets for any rriq ^ 0 


(4.1) 

(4.2) 


2m,a‘Y,(P^(x)P^m „ =-(S”(0)> 


(4.3) 


Although, as the WTI (14.dj) itself shows, there is no mixing between 5° and 
the identity operator with a cubically (or a linearly) divergent coefficient, 
there is still room for a quadratically divergent term proportional to niq. 
Indeed, the l.h.s. of eq. (Q is equal to the piece where intermediate states 
are inserted plus a divergent contribution of the kind mq/a? coming from 
the (integrated) short-distance singularity of the correlator (P^(x)P^(O)) at 
X = 0. This term should be brought to the r.h.s. of the equation, thus leaving 
finite (subtracted) expressions in both sides of the resulting equation. 

We can now repeat on the lattice the argument that in the continuum 
leads to the classical GMOR relation, if we assume that spontaneous chiral 
symmetry breaking occurs in the limiting continuum theory, i.e. if we assume 


that 


S 


lim (G|5°|G) 

mn—»-0 


cont 

mq 


^0. 


(4.4) 


As in the formal continuum theory, we insert a complete set of states in the 
subtracted correlator in the l.h.s. of eq. (I13D. We obtain in this way 


m 


L 

rriq 


»|-(n|s.»„jfi)] 



(4.5) 


where we have explicitely written down only the contribution coming from 
the pion pole. Dots are terms due to the intermediate states that stay mas¬ 
sive as rUq —>■ 0 as well as terms vanishing with rriq faster than linearly (the 
latter include in particular O(a^m^) terms stemming from cutoff effects in 
the sum over xq A- The subscript “sub” is to remind us that it is the prop¬ 
erly subtracted chiral condensate (or, more precisely, the vacuum expectation 

0(a) terms can arise from the sum over ccq in the l.h.s. of eq. H4.,SI) . as in Mtm-LQCD 
symmetry arguments ensure that only O(a^) discretization errors can affect expectation 
values of parity even operators. 
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value of the properly subtracted iso-singlet scalar density operator) that en¬ 
ters this equation. Notice that, as expected, the r.h.s. of eq. (ji3D is a hnite 
renormalization group invariant quantity in the limit a —0. 

Once the leading “IR divergent” cutoff effects have been canceled out, 
the use of the Symanzik expansion in the r.h.s. of eq. (Ii3D yields the formula 


ml± 


= 2m„ 


|(0|P±|7r±)p 


cont 


mq=0 


l + a^J2b,A']+..., (4.6) 

/’\n 


e>o 


where dots denote less dangerous “IR divergent” lattice artifacts compared to 
those explicitely shown as well as contributions of higher order in niq. In get¬ 
ting eq. (BSD we have used the fact that the continuum limit of — 
at vanishing quark mass is S 7^ 0 (see eq. (I4.4j) h 

From the above analysis it follows that, in the region where the series 
in eq. (BH) converges (i.e. at least where the order of magnitude inequality 
ml ~ mq > a? is satished), the squared mass of the charged lattice pion is 
linear in (up to small O(a^) and 0(771^) effects). Thus, vice-versa, we can 
imagine to use deviations from the established linear behaviour possibly seen 
at small mq as a workable criterion to determine the minimal value of mq 
at which simulations can be performed before being set-off by discretization 
effects. 


5 Hadron masses and pion decay constant 

We wish to discuss in this section some peculiar issues concerning the mag¬ 
nitude of the O(a^) discretization artifacts affecting hadronic energies (in 
particular masses) and the pion decay constant. 

5.1 Hadron energies 

As a consequence of the special form of the diagrams contributing to energies 
(self-energy diagrams), one can get convinced that the latter are less “IR 
divergent” than the correlators from which they can be extracted, in the 
sense that at hxed order in a, the most “IR divergent” lattice corrections 
to their continuum limit contain one overall factor 1/m^ less than the “IR 
divergent” cutoff effects appearing in correlators. The reason is that at least 
one among the inserted / £odd operators gets absorbed in a multi-particle 
matrix element, with the consequence that it is not anymore available for 
producing a pion pole 
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An explicit calculation shows that for the difference between lattice and 
continuum energy of the hadron 



one gets at order the estimate 





(ff|4|7r°(0))(7r°(0)a„(q)|4|a„(q)) 

2^an(q) 


cont 


conn 


1 + 0(ml) 


, (5.2) 


where the subscript conn denotes the completely connected part of the matrix 
element (7r°(0)an(q)|f'5|an(q)). In eq. ()5.2j) the “IR divergent” 0(a^/m^) 
piece comes from the continuum correlator with two £odd = nf's + O(a^) in¬ 
sertions. The latter gives also rise to O(a^) “IR hnite” corrections. A further 
O(a^) “IR hnite” correction comes from a single insertion of £even = + 

0(0"^) and contributes a term proportional to a^(Q;„(q)|f'6|Q(„(q))/2T^o„(q)- 

It should be noted that the “IR divergent” lattice artifact in AEa„\a '2 is 
reduced to an “IR hnite” correction after anyone of the two “cures” described 
in Sect. 01 

Specializing the formula (EH to the case of pions, one obtains the inter¬ 
esting result that the diherence between charged and neutral pion (square) 
masses is a hnite O(a^) quantity even if the critical mass has not been set to 
its “optimal” value or the clover term has not been introduced. The reason 
is that the leading “IR divergent” contributions shown in eq. (EH are equal 
for all pions (as one can prove by standard soft pion theorems ng), hence 
cancel in the (square) mass diherence. This conclusion is in agreement with 
results from yPT [16] . 

5.2 Pion decay constant 

Data on a (quenched) computation of /jr, carried out in Mtm-LQCD us¬ 
ing the value of the critical mass obtained from the vanishing of the pion 
mass, show cutoh ehects that are seen to deviate at small bare quark masses 
from the straight line extrapolation drawn from large masses jl]. This be¬ 
haviour (called “bending phenomenon” in ref. |1]) sets in at values of rUg 
around aA^cD> ^QCD rsj 200 MeV. Furthermore the detailed recent scaling 
test of m indicates that the “bending phenomenon” is an O(a^) deviation 
with a magnitude which increases as m^± is lowered. 
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In the works of refs. H Ca ESI El the lattice pion decay constant, 
was extracted from the formula [ 2 ] 


= 2m, 


(f2|P±|7r^) 


9 


(5.3) 


from which it is seen that is a ratio of two lattice quantities. 

It should be noted that in general the matrix element (niP^lTr^^) in the 
numerator is affected by leading “IR divergent” cutoff effects which however, 
like in the case of hadron masses, are multiplied by an extra factor of m^ ~ 
rriq. This peculiar property can be traced back to the invariance of the 
correlator {P^{x)P'^{y))\^^ (from which the matrix element (r2|P^|7r^)|^^ 
is extracted) under axial rotations around the third iso-spin direction. At the 
price of shifting the bare quark mass from rrig to m^ = rrig -|- 0{a?/rrig) one 
can, in fact, always think of having performed the axial rotation with an angle 
such to bring the critical mass to its optimal value. A Symanzik expansion of 
(P^(x)P^(?/))|~ evaluated with the optimal critical mass then shows that 
the leading “IR divergent” cutoff effects (which from this viewpoint arise 
owing to rrig ^ nig only) are softened by the extra multiplicative m^ factor. 
Similarly it follows from the discussion of Sect. 15.11 (which is consistent with 
the previous argument) that also is affected by leading “IR divergent” 

cutoff artifacts softened by an overall m^ factor. The latter, however, drops 
when the relative lattice correction is considered. 

Inserting in eq. (lESD the behaviour of numerator and denominator we 
have just discussed, we hnd that is in general affected by “IR divergent” 
relative corrections that are fully leading in our nomenclature, i.e. of the kind 
(a/m^)^^, k >1. Of course, if the appropriate clover term or the “optimal” 
critical mass is introduced in the action, these “IR divergent” relative lattice 
artifacts are reduced to a? {a? / , k>l effects. 

To be concrete at order a? one finds that in (see eq. (lESD) relative cor¬ 
rections of the type a^/m^ are present, which come only from /rnq)\^^. 
Less important corrections of the type a? /rrig and a? come instead from both 
(r2|P^|7r"f') 1^^ and After the clover improvement or the “opti¬ 

mal” critical mass cure, the relative cutoff effects in get reduced to plain 
a? terms. Results consistent with these features have also been derived in 
yPT studies of Mtm-LQCD jHl EHl at (next-to) leading order. 

A beautiful confirmation of the validity of the analysis presented in this 
section comes from the fact that, when the critical mass is set at its “optimal” 
value (yielding ^°^^{niq) = 0(am^)), no “bending phenomenon” is anymore 
visible in the data, as demonstrated by the results of ref. [EIEIEI (see 
also j2Ij). 
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6 Conclusions 


We have shown in this paper that lattice correlators in Mtm-LQCD are af¬ 
fected by discretization artifacts that tend to become large when the quark 
mass is decreased. Cutoff effects of this kind can be appropriately described 
in terms of the Symanzik LEEA of the lattice theory, and turn out to arise 
from multiple pion poles associated with the insertions of the parity odd (and 
iso-spin non-invariant) piece of the effective action, Jd^y CodA^y)■ 

At order the leading “IR divergent” cutoff effects (those with h = 2km 
eq. (inD) are associated to 2fc-fold pion poles. Since their residues contain 2k 

factors of the matrix element C = |(^|'^odd|7r°(0))| = |a(r2|f'5|7r°(0))|-|-O(a^), 

all these dangerous contributions are eliminated from lattice correlators if we 
can let vanish sufficiently fast as the pion mass goes to zero. This can 
indeed be achieved if Mtm-LQCD is improved a la Symanzik by the inclusion 
of the lattice clover term in the action or by appropriately tuning the critical 
mass parameter to what we have called its “optimal” value. 

In both cases the left-over “IR divergent” cutoff effects are at worse of 
order a^(a^/m^)^“^, k >1. 

The discussion we presented in this work does not make use of chiral 
perturbation theory and treats rUg as a fixed quantity in the limit in which 
a —0. By using the notion of Symanzik LEEA, one is able to obtain results 
to all orders in a. It must be recalled, however, that, as recently pointed 
out in ref. jH], a xPT analysis of (M)tm-LQCD shows that it is possible to 
reabsorb (at any given order in the chiral power counting of choice) a whole 
tower of cutoff effects in a shift of the chiral lattice vacuum, thus extending 
the radius of convergence of the associated expansion 

Our detailed findings are in agreement with the results of refs. BBU 
concerning the quark mass dependence of lattice artifacts and the existence 
of an “optimal” choice for the critical mass yielding a complete cancellation 
of the leading “IR divergent” cutoff effects. In particular, we agree with 
the authors of ref. [S] that, once the “optimal” value of the critical mass 
is employed in Mtm-LQCD simulations (or the clover term is introduced), 
the tight (order of magnitude) limitation rUg > a on the quark mass is sub¬ 
stantially weakened and brought down to the much more favourable bound 
rrig > a?, beyond which large mass-dependent cutoff effects may show up in 
simulations at hxed a. 

We conclude this paper with a remark which explains how the optimal 
critical mass prescription can be viewed as a peculiar mapping between bare 
and renormalized parameters and may be useful in practice to correct for a 

®We are indebted to S. Sharpe and M. Golterman for valuable discussions on this point. 
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not optimal choice of the critical mass. Indeed, the fact that it is possible 
to get rid of all the leading “IR divergent” cntoff effects by working at an 
optimal valne of the critical mass implies that, even if simnlations are carried 
ont with a non-optimal choice of Mcr, one can still effectively eliminate all 
the dangerous {a/iriqY^ [k > 1) terms by means of a slightly more elaborate 
analysis of the available simulation data. This analysis consists in using an 
appropriately redehned expression of the renormalized quark mass (by O(a^) 
terms), accompanied by a small (0(a)) chiral rotation in the third iso-spin 
direction of all quantities that have non-trivial chiral transformation prop¬ 
erties By construction the lattice correlators that are obtained in this 
way will have a Symanzik expansion where the leading ” IR-divergent” cutoff 
effects are just absent, while automatic 0(a) improvement is obviously pre¬ 
served (the deviations from the initial simulation values of the critical mass 
and twist angle are in fact 0(a) quantities). However this analysis requires 
some knowledge (or a simultaneous determination, which can be obtained 
by working at several values of (Mo,mg)) of scale-independent combinations 
of the renormalization constants of operators belonging to the same chiral 
multiplet, such as Zp/Zgo = ZpZ^ or ZyjZ^- 

Indeed the prescription we have described is the microscopic counter-part 
of what is suggested in ref. jH] at the level of the effective chiral theory in 
order to reabsorb the contributions of these lattice artifacts via a proper 
redehnition of the vacuum state. 

Being based on the assumption that an acceptable (i.e. 0(a) accurate) es¬ 
timate of the critical mass is already known, this remark is of no use to obtain 
such an estimate (the argument would become circuitous). It should also be 
observed that the effort of determining by how much the renormalized quark 
mass have to be shifted and chiral operators rotated is essentially equivalent 
(modulo the knowledge of Zp/Zgo and analogous renormalization constant 
ratios) to enforce the optimality condition in the way we have described it 
in sect. Id.21 
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®This is the kind of procedure we have alluded to in sect. E2 below eq. (lot . Detailed 
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relation between bare and renormalized mass parameters and currents (or more general 
operators) in tmLQCD. 
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Appendix - improvement of Mtm-LQCD 


We prove in this Appendix the absence of all the terms of order ^ > 0, 

in the Symanzik expansion of the expectation valne of parity-even (mnlti- 
local) gauge-invariant and multiplicatively renormalizable operators in Mtm- 
LQCD. Here to be precise by “parity-even” we mean an operator which 
goes exactly into itself under parity. This result follows from the symmetry 
P X Vd X {fT^q —iTT'q), eujoyed by the lattice action, where {xp = (—x, t)) 


Uo{x) 

->■ Uq{xp) , 

Uk{x) 

- Ul{xp- 

^_{x) 

7o^(xp) 

^ i/>(x) 

'0(xp)7o 


is the physical parity of the theory, and 


Vd-. 


Uf,{x) Uj^i-x-afi), 
-ipix) —*• x) 

ip{x) —»• x) 


A: = 1,2,3 


/i = 0,1,2,3 


(A.l) 


(A.2) 


The proof of improvement of Mtm-LQCD follows immediately 

from the observation that necessarily the Symanzik LEEA which describes 
the lattice artifacts of Mtm-LQCD is invariant under (the continuum version 
of) the transformation P x Vd x {irtq —> —niq). This invariance implies, in 
particular, that all the terms of order /c > 0, in (|22D are odd under 

parity. In fact, Vd x {niq —>■ —mq) counts the parity of the overall dimension, 
d, of any product of elementary fields and mass factors, by multiplying it by 
the phase factor exp{ind) (besides inverting its space-time argument). As 
a result, all the continuum correlators in the Symanzik expansion that are 
multiplied by an odd power of a necessarily correspond to expectation values 
of parity odd operators. Since parity is an exact continuum symmetry, they 
all vanish. 

Few comments are in order here. 

1) It is interesting to observe that automatic ©(a^^"*"^). A: > 0, improve¬ 
ment (in the Symanzik sense) is a very robust property of Mtm-LQCD. Being 
based on straight symmetry arguments, it is independent of any considera¬ 
tion on possible phase transition scenarios piin]. In fact, for any value of 
niq, sufficiently close to the continuum limit the chiral phase of the vacuum 
will be finally determined by the quark mass term in the action and not by 
the Wilson term. 
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2) At maximal twist the r-parity properties of the critical mass are imma¬ 
terial to the above argument for automatic fc > 0, improvement. 

The reason is that in this situation (but not for generic values of the twist an¬ 
gle) the Wilson term, and thus necessarily also its critical mass counter-term, 
are odd under parity. This clearly remains true for whatever determination 
of the critical mass is taken, independently of its behaviour under r —>■ —r. 

On the contrary, at non-maximal twist, the proof of 0(a) improvement 
(via Wilson averaging) requires the critical mass to be an odd function of r. 
As explained in the papers of ref. [2j, this condition can (and should) always 
be imposed when determining Mcr(r). 

3) Finally we note that the conclusions about automatic 0(a^^’''^) im¬ 
provement reached in this Appendix can be extended to the case of maxi¬ 
mally twisted quarks with non-degenerate masses jS] . If the physical basis of 
ref. |S] is employed, a proof completely analogous to that given above for mass 
degenerate quarks can be given relying now on the symmetry P x Pd x {fUq —> 
—mq) X (cg —>• —Cg), with the bare mass splitting within the quark pair. 
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